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INTRODUCTION 



I . 

Search and Rescue (SAR) is one of the United States 
Coast Guard's most visible and well known primary missions. 
The backbone of the Coast Guard's SAR effort is the network 
of small boat Stations distributed at various points along 
our nation's shoreline so as to provide adequate coverage 
for coastal boating. The bulk of these Stations’ SAR fleets 
consists of two boat types, the forty one foot utility boat 
(UTB) and the forty four foot motor surf boat (MSB) . These 
two hull types account for virtually all the short range 
surface SAR response provided by the Coast Guard. 

While it would be ideal in terms of public safety to 
maintain an unlimited number of SAR resources available for 
immediate response, such a policy is economically infeasi- 
ble. This gives rise to the obvious though difficult prob- 
lem of determining the optimal allocation of a scarce 
resource, in this case, boat crews. The allocation of boat 
crews vice boats is discussed since it is generally Coast 
Guard policy to allow at least one more boat than boat crew 
per Station to increase the likelihood of having both a boat 
and crew available for response simultaneously. Thus, the 
boat crew is the limiting resource. 

The Search and Rescue Division of the Office of Opera- 
tions, Coast Guard Headquarters, maintains ultimate responsi- 
bility for this task. In deciding on each Station's manning 
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level. Division personnel are interested in the effects 
felt by local SAR systems caused by varying available re- 
sources. Specific effects can be measured by such measures 
of effectiveness (MOE's) as the probability of a delayed 
distress response caused by occupied boat crews, the average 
duration of the delay, and the average number of distressed 
craft having to wait. Obtaining such information would 
surely assist planners in arriving at some satisfactory 
allocation . 

These MOE's suggest the use of queueing techniques to 
analyze the operation of SAR resources at a Station. The 
incidence of distress cases can be considered as customer 
arrivals. Upon entering the system, cases may either wait 
for service or receive service immediately upon finding an 
idle boat crew. The length of time a boat is deployed on a 
given case will comprise the service time. Information on 
both arrival times and service times is contained in the SAR 
data base maintained at Coast Guard Headquarters. The data 
base consists of information recorded by each unit on the 
Search and Rescue Assistance Report (CG-5151) [Ref. 1] each 
time they respond to a case. 

This thesis is concerned with examining the available 
data in order to determine an appropriate queueing system 
which will best model the situation, and then providing solu- 
tions for the system in terms of some MOE's which could be 
of value to the SAR planner. Chapter II of this thesis 
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contains a description of the manner in which distress cases 
arrive at a Station. Justification for considering only a 
finite time window in order to identify a period of homo- 
geneous arrivals during the busiest SAR season is presented, 
along with the results of the data-fitting analysis. A 
discussion of service times follows in Chapter III, again 
along with the results of efforts to find an appropriate 
distribution to fit the data. The samples were taken from 
four environmentally diverse Stations intended to represent 
a reasonable cross section of Coast Guard operating areas. 

The information was drawn from the fiscal year 1980 data base 
as it was the most recent available at the time of commencing 
this study. Chapter IV provides a list of desirable system 
descriptors of use to a Coast Guard SAR planner and equates 
these to some common measures of effectiveness often used 
in queueing analysis. Two practical computational methods 
for calculating the desired information and some work con- 
sidering the introduction of a two class priority scheme are 
then presented in Chapter V. The final chapter contains a 
brief summary of the work. 
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II. DISCUSSION OF THE ARRIVAL PROCESS 



A. DEFINITION OF AN ARRIVAL 

The definition of an "arrival" warrants some discussion. 
From the time a vessel perceives its distress to the moment 
a Coast Guard vessel arrives alongside, several communica- 
tions have transpired, both between the distressed vessel 
and the Coast Guard, and within the Coast Guard's SAR struc- 
ture. Exactly when during this time a case is said to 
"arrive" needs to be decided. 

A distressed craft will communicate its situation by 
any means available. This could be by a mayday message 
broadcasted over a radio if so equipped, by one of many 
visual distress signals such as lighting a fire on deck or 
sending up a red flare, or maybe simply by shouting to 
another vessel in close proximity. This initial call for 
help may or may not fall directly on Coast Guard ears, but, 
in most cases, the Coast Guard's first notification will 
arrive at a Radio Station. While not normally equipped with 
life saving resources, they do act as a communications 
facility and can effectively relay relevant information to 
the cognizant Rescue Coordination Center (RCC) . The RCC 
acts as a command center. Here the decision is made as to 
the appropriate resource to dedicate to the mission. If it 
is felt the best response could be provided by a Station's 
UTB or MSB, that Station is contacted, either by radio-teletype 
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or telephone, and ordered to dispatch on the case. The 
Station then assigns a boat and a crew and gets them under- 
way towards the vessel in distress. 

Since only those data recorded on the SAR Assistance 
Report [Ref. 1] are available, three entries emerge as 
candidates; day/tirae CG notified, day/time unit notified, 
and day/time underway. Of the three, the second, day/time 
unit notified, appears to be the most appropriate. The 
purpose of this work is to gain insight into the conse- 
quences of varying resources at the Station, other factors 
remaining constant. It therefore would not be appropriate 
to consider a case as having arrived prior to the Station 
being made aware of its need for response, eliminating the 
first candidate. As for the third, day/time underway would 
exclude time spent in the queue awaiting an idle resource, 
while it is an aim of this thesis to use that quantity as 
an MOE. This leaves day/time unit notified as the logical 
time of arrival. 

B . DATA ANALYS IS 

1. Homogeneous Time Window 

As suggested in the Introduction, Station SAR demand 
is far from homogeneous, due largely to the fact that the 
"customers" serviced are primarily members of the recreational 
boating population. Experience and the data clearly show 
SAR case load and boating density are lightest during the 
colder, more severe winter months and during the normal five 
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day work week. Rather than consider these non-homogeneous 
arrivals in the queueing analysis, an attempt was made to 
find a finite time window during the busiest hours of the 
busiest days of the busiest month in which arrivals could 
be considered to be occurring at a constant rate. 

2 . Statistical Procedure 

Graphical displays of weekend SAR data suggested the 
period 0900 to 1800 on Saturday and Sunday as a candidate 
for a time interval over which cases arrive at a constant 
rate. Times from the beginning of the time window to the 
first case arrival and from the last arrival to the end of 
the window were not considered. Efforts to fit these inter- 
arrival times, drawn from four environmentally diverse 
Stations (New Haven, Ct., Rockaway, N.Y., Cape May, N.J., 
and Coos Bay, Ore.) , to the Exponential distribution proved 
successful. July data was chosen to represent the season 
of busiest SAR activity. The decision to accept the Exponen- 
tial distribution was based on Chi-square goodness of fit 
tests at the 95% level, using the first moment estimator, 
the reciprocal of the sample mean interarrival time, as the 
estimate of the arrival rate A. Table I contains a summary 
of the results of the tests on each Station's interarrival 
data. The computed Chi-square statistics are all well below 
the corresponding critical values, providing no evidence to 
reject the assumption that the data came from an Exponential 
population. The acceptance, therefore, of the Exponential 
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TABLE I 



Exponential Chi-Square Test Results 



Station 


A (hrs.)"^ 


^computed 


^ . 9 5 , n- 1 


Cape May 


1.797 


3.30 


12.59 


Coos Bay 


1.022 


2.37 


12.59 


New Haven 


1.018 


0.35 


9.49 


Rock aw ay 


1.340 


3.22 


14.07 


distribution 


implies that the 


Poisson Process is 


an appro- 



priate descriptor of the arrival phenomenon. 
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III. DISCUSSION OF SERVICE TIMES 



A. DEFINITION OF SERVICE TIMiE 

As with arrivals, exactly what is to be considered 
service time needs to be defined. Unlike a bank or barber 
shop, customers do not physically arrive at a well-defined 
facility and commence being administered to, face to face, 
by the first available server, nor do customers "queue up" 
in a waiting line when awaiting service. In most instances, 
a typical SAR case occurs some distance offshore. The crew 
must get a normally moored vessel underway, transit the 
distance from the Station to the distressed craft, render 
necessary assistance, and, in more cases than not, tow the 
disabled vessel to safety. Since all elements of the afore- 
mentioned procedure are part of servicing a case, all should 
contribute to service time. 

A considerable problem arises due to the nature of the 
available information. Data recorded include day/time unit 
notified and day/time underway as previously discussed, as 
well as sortie time. This latter term is defined in the 
Coast Guard piablication Search and Rescue Data System Manual 
(COMDTINST M5230.10) [Ref. 2] and basically involves the 
underway time of the rescue vessel. More specifically, 
sortie time is either the time from departure until the time 
of return to homeport, or, if diverted from another mission, 
the time of diversion to the time of return to the original 
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mission. The problem arises in considering the elapsed time 
from notification to time underway. Given an available 
resource (boat) , this time typically involves either a short 
walk from an adjacent shore facility to the boat or often a 
radio call directly to the boat if the crew is already 
aboard. However, if a case were queued, the waiting time 
would accumulate here with no accompanying indication that a 
resource was unavailable. In other words, there is no way to 
isolate that portion of the time from notification to the 
time underway due to normal startup time, which of course 
should be included in service time, from that time spent 
waiting for the availability of a resource. Examination of 
the raw data, however, revealed that the time from unit 
notification to time underway, which would include any 
queueing time, was no more than ten minutes 90% of the time. 
In further analysis, therefore, sortie time will be used as 
a surrogate of service time. 

B. DATA ANALYSIS 

1 . Parameter Estimation 

Sortie times for the four Stations were then dis- 
played in stemleaf plots to aid in discovering an appropriate 
distribution to fit the data. Only those service periods 
begun during the homogeneous time window were considered. 
Examination of the displays consistently suggested the use 
of the Gcimma or Erlang distributions, of which the latter 
was chosen due to the availability of many queueing results 
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involving its use. Due to the presence of two parameters in 
the Erlang distribution, however, three methods of parameter 
estimation were investigated. Letting Y be a random variable 
representing service time, the appropriate PDF is 

,1,K 1 

f(y) = - ( ■ ]<- 1) ! y ® ^ ' y > 0 (3.1) 

2 

with mean kg and variance kB . 

The first alternative involved using the method of 
moments, setting the sample moments equal to the correspond- 
ing distributional moments; 

Y = kS and = kB^. 

This yielded the parameter estimates: 




Since k must take on an integer value, k must be rounded 
to the nearest integer . 

The second method varies from the first only in its 
estimate for B and accounts for the rounding error created 
when rounding k to the nearest integer. The estimate for 
B is therefore recalculated using the integer valued estimate 
of k and the relationship 
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mean 



kS 



yielding 



k 

The third method investigated is presented by Read 
[Ref. 3] and is derived through the Gamma distribution, then 
making the appropriate conversions to the Erlang distribu- 
tion. His estimates for the Gamma parameters 

a = and 6 = Y - - 

Y2-1 ^ 



where 



z 



1 ^ 
- y 

n . 



i=l 




are said to be asymptotically more efficient than the simple 
method of moments estimators for k > 3, 

2 . Statistical Procedure 

The Erlang parameter estimates for each of the four 
Station samples were calculated using each of the aforemen- 
tioned methods. Using the estimates provided. Chi-square 
goodness of fit tests were conducted. The second method 
consistently produced the smallest Chi-square statistics, 
well within a 95% acceptance region, and thus was chosen to 
provide future parametric estimates. The results are 
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displayed in Table II. It should be noted that the perfor- 
mance of the third method in terms of the Chi-square tests 

/V 

was not surprising due to the small values of k generated. 
Perhaps in future application, if samples produce larger 
estimates of k, the use of this method would prove to be a 
more competitive alternative and should therefore not be 
completely abandoned. 



TABLE II 

Erlang Chi-Square Test Results 







Method 


#1 




Method 


#2 




Method 


#3 




Station 


k 


6(hrs.) 


2 

Y 

'^oonp 


k 


B (hrs . ) 


2 

\xfftp 


k 


6 (hrs . ) 


2 

^oonp 


2 

^.95,n-l 


Cape May 


2 


.703 


10.16 


2 


.778 


8.24 


2 


.668 


21.01 


12.59 


Coos Bay 


1 


1.169 


4.27 


1 


1.051 


2.18 


2 


.573 


11.24 


9.49 


New Haven 


2 


.810 


U.84 


2 


.689 


6.70 


3 


.534 


11.16 


9.49 


Bockaway 


2 


.620 


7.39 


2 


.577 


6.14 


2 


.611 


5.39 


11.07 
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IV. MEASURES OF EFFECTIVENESS 



The results from the two previous chapters support the 
contention that the M/E^/c queue may be of value in des- 
cribing the operation and allocation of Search and Rescue 
resources at a Coast Guard Station. It remains to determine 
a method of extracting needed measures of effectiveness to 
aid in evaluating the relative merits of competing resource 
allocations . 

In administering the Search and Rescue program, Coast 
Guard planners must decide upon, as previously discussed, the 
number of ready boat crews to allocate to each Station. 

Their problem is that of supplying the customers, in this 
case the recreational boating public, with some satisfactory 
level of SAR capability while attempting to minimize system 
costs. Certain relevant measures of system performance 
obtainable through queueing techniques could prove of great 
utility in this endeavor. These measures are described 
below . 

A. PROBABILITY OF A DISTRESS CASE WAITING FOR SAR SERVICE 

An important yardstick with which to measure the adequacy 
of a given resource level is the estimate of the probability 
that a vessel in distress will have to wait for an idle res- 
cue craft. Due to the emergency nature of the required 
response, it is of paramount importance to maintain a con- 
sistently high readiness posture. Unlike a barber shop or 
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gas station, the consequences of waiting for service may be 
grave and even potentially fatal . A method through which 
one can assess, for varying boat crew allocations, the likeli 
hood of failure to meet the SAR demand at once during a 
peak period would be of great value. Using queueing analy- 
sis, the quantity termed delay probability provides this 
result and will be calculated in the next chapter. 

B. AVERAGE NUMBER OF DISTRESS CASES WAITING FOR SERVICE 

AND IN THE SYSTEM 

The analyst might choose to weigh the relative merits of 
two or more proposed allocation schemes based on the number 
of distressed vessels, on the average, having to wait for an 
idle Coast Guard crew, or similarly, on the average number 
of active cases, both those being assisted and those waiting 
for service, at any time. Clearly, the ordinal ranking of 
fewer is better would be appropriate. The queueing expres- 
sions L (Average queue length) and L (average number in 
system) provide methods of calculating these averages. 

C. AVERAGE TIME SPENT WAITING BY A DISTRESS CASE FOR AN 

IDLE CREW AND AVERAGE TIME IN SYSTEM 

The planner might wish to provide boat crews sufficient 
to keep the average duration of a distressed case's delay 
below some maximum value. The total time, on the average, a 
case remains active, including both waiting and service time, 
might provide another MOE of similar value. These quantities 
designated W^ (average waiting time) and W (average time in 
system) are readily obtainable through queueing techniques. 



21 



D. DISTRIBUTION OF WAITING TIMES OF A DISTRESS CASE 

Often a SAR system manager would be interested in the 
likelihood that a distressed vessel's delay exceeds some 
fixed length of time. This type of analysis lends itself to 
the establishment of some threshold duration or maximum 
acceptable delay with very high probability. Queueing 
analysis could allow the planner to determine, for instance, 
the number of boat crews required to ensure that the likeli- 
hood a vessel must wait more than a half hour for an idle 
rescue boat does not exceed ten percent. A knowledge of the 
distribution of waiting times would therefore serve as 
excellent input to the decision making process. 



22 



V. COMPUTATIONAL METHODS 



This chapter presents two reasonably simple methods for 
obtaining the aforementioned MOE's given the appropriate- 
ness of the M/Ej^/c queue as a modelling vehicle for the 
Station SAR operation, and then extends the results to 
consider a priority scheme other than FIFO. 

The work thus far has concentrated on providing neces- 
sary background on the Coast Guard SAR system, finding 
characterizations of the arrival process and service times, 
and selecting MOE's. Prior to presenting the two computa- 
tional methods previously discussed, a more detailed descrip- 
tion of the M/Ej^/c queueing system is in order. 

A. DESCRIPTION OF THE M/Ej^/C QUEUE 

The M/Ej^/c queue is characterized by a homogeneous 
Poisson arrival stream with rate A whose interarrival times 
have PDF 



f (x) 



A e 



-Ax 



X > 0 



and Erlang service times with PDF 



f (y) 
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The FIFO (first in, first out) priority discipline is 
assumed and the following variable and parameter definitions 
apply. 

c = total # of servers 

X = arrival rate = reciprocal of mean inter- 

arrival time 

k = Erlang shape parameter 

6 = Erlang scale parameter 

U = service rate = reciprocal of mean service 
time 

p = X/c\x = system intensity 

N = total number in system 

N = total number in queue 

q 

L = E[N] = mean number in system 

L = E [N ] = mean number in queue 

q q 

T = total time in system 

T = total time in queue 

q 

W = E[T] = mean time in system 

W = E [T ] = mean time in aueue 

q q 

B. TABULAR METHOD 
1 . Background 

The theory surrounding the M/E, /c queue has been 
sufficiently developed to a point where several numerical 
methods have been shown to produce adequate steady state 
results. Due to the complexity of the computational require- 
ments, however, direct utilization by a layman in a decision 
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making capacity (in this case, a Lieutenant in the Search 
and Rescue Branch) would be impractical, if not impossible. 
Hillier and Lo [Ref. 4], in an effort to make complex results 
readily obtainable by would-be users, have prepared, for a 
finite though practical subset of possible values of system 
parameters k, c, and p, a series of tables presenting many 
desirable measures of steady state system performance. They 
cite recent computational treatments by Mayhugh and McCormick 
[Ref. 6], Mayhugh [Ref. 7], and Heffer [Ref. 8], and the 
comprehensive analysis of the E^/Ej^/c system by Yu [Ref. 9] 
as the basis for their tables. Of specific interest for 
this application is their Table I, partially reproduced 
herein as Appendix A, which provides the desired measures of 
effectiveness for the following parametric combinations. 

k = 2 ; c = 2 , 3 , . . . , 10 

k = 3; c = 2, 3, 4, 5 

k = 4; c = 2,3 

k = 5, 6, 7, 8; k = 2 

p = .10 , .20, . . . , .50 , .55, . . . , .95, .98, .99 

2 . Use of Tables 

Use of the tables requires entering arguments c, k, 
and p. The parameter c, number of boat crews, will most 
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likely be the design variable whose determination is the 
object of analysis. Estimates for the parameters k and 
p will be provided by the SAR data from the Station under 
study. Recommended parametric estimates, as previously 
discussed, are 

A = 1/X 

'' —2 2 

k = y /S then rounding to nearest integer 

/N /\ 

3 = y/k 

/N ^ /S 

P = 1/kS 

/\ /N 

p = A/cp 

The reader is referred to Chapters II and III for discussions 
concerning the applicability of the SAR data base and the 
criteria for the determination of the homogeneous time window 
to satisfy the assumptions of the M/Ej^/c queue. 

Since results are listed for discrete values of the 
system intensity p, interpolation will undoiobtedly be neces- 
sary. Since no interpolation scheme is suggested in Ref. 4, 
and since visual examination of the tabular entries suggests 
it would not be unreasonable, the linear technique is assumed. 

The aforementioned measures of effectiveness are 
obtainable, directly or indirectly, from the tables. The 
probability of requiring a distress case to wait for response 
due to preoccupied rescue boats , otherwise called delay 
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probability, is listed separately at the end of each section 
of parametric combinations. Though interpolation for the 
exact value of system intensity p is required, this impor- 
tant result is otherwise available at the analyst's finger- 
tips. Likewise, average number of cases waiting (L ) and 
average number of cases in the system (L) are similarly 
accessible. Calculation of average time spent waiting for 
service (W^) and average time in the system (W) is made 
possible through the well known Little's [Ref. 10] relation- 
ships 



W = and W = ^ . 

q X X 

In addition, though not directly tabulated, Hillier 
suggests an approximation method for estimating the distri- 
bution of waiting time T . The derivation is based upon the 
observation that an Erlang random variable with mean 1/u may 
be expressed as the sura of k Exponential phases, each with 
mean l/k^, and that, when there are no idle servers, comple- 
tions of service phases occur according to a Poisson Process 
with rate cky . This leads, as discussed in Ref. 4, to the 
following probability distribution function. 

OO 

P[T > t] = y P [N=n]P[D < k(n-c+l) -1] , t > 0 
^ n=c 

where D is a Poisson random variable with mean ckpt. Since 
the first term is tabulated by Hillier and Lo in Ref. 4, and 
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cumulative Poisson probabilities are readily obtainable, 
either in tabular form or through direct computation, the 
entire expression is easily computable. 

3 . Examples 

By manner of illustration, the following two tables 
demonstrate the effects on system behavior resulting from 
varying the number of servers (boat crews) in an otherwise 
stable queue. Table III contains results stemming from the 
author's arbitrary and hypothetical selection of system 
parameters A and y, while Table IV displays results produced 
from actual SAR data drawn from Coast Guard Station Rockaway , 
N.Y. 

TABLE III 
Hypothetical Case 

A = 1.5/hr k = 3 y = 1.0/hr 





c = 2 


c = 3 


c = 5 


MOE 


p = .75 


o 

LO 

ti 

CL 


p = .30 


P [N = 0] 


.139 


.208 


.223 


P[T >0] 

q 


.6 39 


.232 


.020 


L 


2.803 


1.667 


1.507 


L 

q 


1.30 3 


.167 


.007 


W (hrs.) 


1.869 


1.111 


1.005 


W (hrs.) 

q 


.869 


.111 


.005 


P [T > .5 


.513 


.105 


.004 



hrs . ] 
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TABLE IV 



Station Rockaway 

X = 1.340/hr k = 2 6 = .577 hrs y = .867/hr 



MOE 


c = 2 
p = .77 


c = 3 
p = .52 


c = 4 
p = .39 


P [N = 0] 


.128 


.195 


.209 


P [T > 0 ] 

q 


.668 


.256 


.084 


L 


3.275 


1.782 


1.605 


^q 


1.735 


.222 


.045 


W (hrs.) 


2.444 


1.330 


1.197 


W (hrs.) 

q 


1.295 


.165 


.033 


P [T > . 5 hrs . ] 


.561 


.088 


.008 



The above tables illustrate quite graphically the 
effects of varying the number of servers or boat crews and 
could justifiably provide the SAR planner with valuable 
information upon which to base a decision. Given the exis- 
tence of some design criterion specifying, for example, a 
maximum acceptable delay probability, this method of analysis 
would clearly indicate the minimum number of boats required 
to achieve it. For example, if Coast Guard policy required 
that the chance of a delayed response not exceed, say, .30, 
Table IV clearly indicates at least three boat crews would 
be required at Station Rockaway in order to meet that require- 
ment during the peak load periods. Similarly, a system re- 
quirement that the average amount of time spent waiting for 
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service be less than some amount, say one half hour, could 
very simply be satisfied through ensuring that the required 
mamber of boat crews are assigned. 

C. COMPUTATIONAL METHOD 

1 . Background 

In his Master's Thesis, Celayir [Ref. 5] developed a 
relatively simple analytic solution to the M/E-^/c queue, 
which has the advantage of supplying rapid, computable 
results without the necessity for extensive tables or 
sophisticated computational capability. Additionally, the 
method is not constrained to the discrete set of system 
parameters k and c covered in Lo [Ref. 4], and thus provides 
a valuable alternative when existing SAR data produce param- 
eter estimates beyond the values mentioned in the previous 
section. Further, Celayir 's results were made plausible 
through comparisons made with three other sources, specifically, 
simulation output, Hillier and Lo's [Ref. 4] tables, and 
observed similarities with the tractable M/M/c queue. The 
reader is referred to Ref. 5 for a more in-depth discussion 
of the derivation. 

2 . Computational Procedure 

Like the previous method, system parameters k, B, 

U, and p are required for the calculations and should be 
estimated from the SAR data base as discussed in Chapters 
II and III. One additional parameter, offered load, is 
mentioned by Celayir and is defined as the ratio of the 
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arrival rate X to the service rate y and can be expressed 
as the product of the number of crews c and the system 
intensity p yielding 



a 



X 

U 



c p . 



It was observed by Celayir that, for equal values 
of system intensity p, delay probabilities (P[T^>0]) for 



the M/Ej^/c queue were virtually identical to those of the 
M/M/c queue. Thus, the delay probability for the latter 
system provides an excellent approximation for that of the 
former and is computable [Ref. 11] as 



This expression, though in a different form, is algebraically 
equivalent to that presented in Ref. 5. Should this calcu- 
lation prove too cumbersome, as may be the case for c > 4 
(although, in this application, a Station is rarely billeted 
for more than three boat crews) , Appendix B contains two 
charts which provide delay probabilities for the M/M/c 
queue given the offered load a as the entering argument. 

Either of these two methods provides a convenient and accurate 




c(A/y) ^ 



(5.1) 



c! (cy - A) PQ 



where 



PO 
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approximation, without the need for interpolation, for the 
probability a distress case will need to wait for an idle 
boat crew . 

An approximation for the average number of vessels 
waiting for service, L^, is derived in Ref. 5 and is based 
upon extensions of the theoretical work in Lo [Ref. 4] . The 
formula stems from Hillier and Lo's approximation 



L 

q 



|(1 + g) (1 + |) 



L 

^1 



where L is the average queue length for the corresponding 
M/M/c queue and the function g is a multiplicative correction 
factor developed by Lo [Ref. 4] to be applied in the compu- 
tation of when extrapolating the tabular results for 
large values of the parameters k or c. But, for smaller 
values of system parameters k and c, g tends to zero. Noting 
this and the fact that Lq^^ may be simply expressed as a 
function of the previously derived delay probability, that is 



_ (A/y) Xy 

'^1 (c-l)l(cy-A) 



2 Pq 



• P [T > 0] , 
cy “A cj 



the average queue length may then be written as 



L 

q 



i( 1 + i) ( h 

2^-^ ^ k^ ^cy-A 



) *P [T > 



01 . 
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It then follows, again using Little's [Ref. 10] relationship, 
that the average time spent waiting for an available crew 
equals 







P[T > 

q 



0 ] 



The average number L of cases in the system at any time and 
the average time W spent in the system are then expressable 
as 



L = 




X 

y 



and 



W = 




L 

X * 



Thus, given system parameters X, k, y, and c and having 
obtained the delay probability P [T >0], approximations for 

q 

average number waiting, average time spent waiting, average 
number in system, and average time in system are all easily 
computable . 

Finally, Celayir [Ref. 5] proposes a direct method 
of calculating the probability that the duration of a dis- 
tressed vessel's wait exceeds a specified length of time. 
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He suggests that the cumulative distribution function (CDF) 



of waiting times is 



(t) = 1 - P[T > 0] e~^^ , t > 0 

q ^ 



where 



b 



2 c y ( 1 - 

-I 



P) 



The interested reader is referred to his work for an explana- 
tion of the derivation. One should note however, that all 
component parts of the expression are obtainable as previously 
described. If the SAR planner, for example, wished to know 
the probability that a potential case would have to wait 
longer than one hour for an available Coast Guard response, 
he would need only to compute 



P[Tg 



> 0 ] 



-b 

e 



Similarly, should he wish to know the chances that an arriving 
case might wait longer than one hour but less than two, the 
expression 



P [Tg > 0] (e"^ - e"^^) 

would provide the answer. It should be emphasized that, 
while Celayir [Ref. 5] showed that results obtained using 
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this expression compared quite favorably with those from his 
simulation and from Hillier and Lo's tables, considerably 
less effort was required to compute his formulas. 

3. Examples 

The following two examples are presented to provide 
a comparison with Tables III and IV. All relevant system 
parameters were held constant. 



TABLE V 

Hypothetical Case 

X = 1/5/hr. k = 3 y = 1.0/hr. 

c = 2 c=3 c = 5 



MOE 


p = .75 
a = 1.50 


p = . 50 
a = 1.50 


p = .30 
a = 1.5( 


P[T >0] 


.642 


.238 


.019 


q 


L 


2. 784 


1.659 


1.505 


L 


1.284 


.159 


.005 


q 


W (hrs.) 


1.856 


1.106 


1.003 


W (hrs.) 

q 


.856 


.106 


.003 


P [T > . 5 hrs . ] 


.441 


.077 


.001 


q 
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TABLE VI 



Station Rockaway 

X = 1.340/hr. k = 2 S = .577 hrs . \i = . 867/hr. 

c = 2 c=3 c = 4 



MOE 


p = .77 
a = 1.54 


p = .52 
a = 1.56 


p = .39 
a = 1.5 


P [T > 0 ] 

q 


.669 


.254 


.077 


L 


3.252 


1.749 


1.582 




1.706 


.203 


.036 


W (hrs.) 


2.427 


1.305 


1.181 


W (hrs.) 

q 


1.273 


.151 


.027 


P [T > .5 hrs . ] 

q 


.513 


.195 


.019 



As can be seen, Tables V and VI resulting from calcu- 
lations using Celayir's method compare quite favorably with 
the figures in Tables III and IV. Certainly, in the context 
of providing the SAR planner with a guide with which to com- 
pare the relative performance of a Search and Rescue Station 
given various numbers of boat crews, both methods are of 
great utility. 

D. PRIORITY CONSIDERATIONS 

The reader is reminded that the queueing techniques 
applied thus far have assumed the FIFO priority discipline. 
This assumption may carry with it far reaching ramifications, 
for, in practice, resources are assigned to distress cases 
according to a continuous priority scale. Each arriving 
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case is evaluated and assessed as to its severity relative 
to other active cases. If, in the judgement of the cogni- 
zant Rescue Coordination Center, the newly arrived distress 
case provides a greater or more imminent threat to life or 
property than other vessels currently being assisted, ongoing 
assistance will be terminated and that resource redirected 
towards the more urgent situation. A possible scenario fit- 
ting this description might be a Coast Guard ULB diverting 
from retrieving a disabled pleasure craft at anchor in a 
protected waterway in favor of rendering assistance to a 
vessel on fire or in danger of going aground. 

1. Background 

Extensive theoretical work concerning the inclusion 
of priority schemes other than first come first served 
exists, but most is concentrated on the more common M/M/1 
queue, that is single server with Exponential service. 
Extensions to systems involving multiple servers and service 
distributions other than Exponential are rare as well as 
extremely complex, and clearly stand as an area for further 
research. The most thorough treatment encountered by the 
author can be found in Jaiswal [Ref. 12] , and the reader is 
referred there for a proper initiation into the subject. By 
way of an overview, however, the following discussion is 
provided. 

When considering priority schemes other than FIFO, 
arrivals are usually thought to belong to one of several 
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r 




priority classes. Prior to the saturation of idle servers, 
customers are served as they arrive. Once all servers are 
busy, however, customers assume a place in the queue accord- 
ing to their priority class, the customers belonging to the 
higher priority classes moving ahead of those in lower classes. 
It can be seen, then, that the potential exists for customers 
of a low priority class to be bumped indefinitely. The 
nature of the priority discipline may be further refined by 
considering whether priorities are preemptive or non-preemptive . 
In the latter, an arriving customer belonging to say the 
highest priority class will move to the front of the queue 
and wait for the next available server even if customers of 
a lower class are currently in service. If, however, the 
priority scheme is said to be preemptive, current service on 
a lower priority customer will be terminated upon the arrival 
of a customer belonging to a higher priority class. In the 
context of a Coast Guard SAR Station, distress cases would 
be considered as having preemptive priorities with the capa- 
bility of interrupting ongoing assistance to a case belonging 
to a lower class . 

2 . Two Priority Class Case 

It has been suggested by personnel from the Search 
and Rescue Branch that two priority classes, labeled serious 
and routine, would suffice. This designation arises from 
information regarding the severity of a case recorded, after 
the fact, on Ref. 1, and therefore contained in the SAR data 
base. Thus, separate arrival rates and candidate descriptive 
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distributions for both classes are readily obtainable. While 
this type of analysis is a gross approximation to the con- 
tinuous set of actual priority classes, it would probably 
provide an improvement over considering the FIFO discipline 
alone. 

While the complexities introduced by priority con- 
siderations create a considerable barrier in obtaining 
exact solutions to the M/Ej^/c queue, one important measure 
of effectiveness, the delay probability, is somewhat tracta- 
ble. If one were to restrict the number of priority classes 
to two, for instance, estimates for the delay probability 
for both are derived for the cases c = 1 and c = 2. 

Considering first the case of the higher priority or 
serious class of distress cases, the direct application of 
techniques discussed in Chapters II, III, and V is in order. 
This stems from the observation that, given the preemptive 
priority discipline, a server occupied with a routine case 
is identical, from the point of view of an arriving serious 
case, to an idle server. Thus, the analyst could determine 
the delay probability of a serious case by considering only 
the arrivals and service times of cases belonging to that 
class from the initial data base, and then proceeding as 
before, while ignoring all routine cases. This would involve 
fitting this abbreviated data to arrival and service distri- 
butions, and then calculating the desired quantities from 
the resultant queueing system. Although unable to confirm 
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it (available data contains no more than six serious arrivals 
per Station ) , the author has seen no evidence to suggest 
that the M/E-^/c system would not again be appropriate. 

The question of the delay probability for the non- 
serious or routine case presents a somewhat more complicated 
problem arising from the fact that a routine case might be 
delayed even after its service has begun. Specifically, if 
all servers are busy, an arriving serious case could preempt 
service on a routine case with the latter assuming a place 
in the queue at the front of other cases belonging to the 
same class. The derivations that follow require the intro- 
duction of some new notation. 

D = total delay time for a routine case; 

= initial delay encountered by a routine case 
upon arrival before its service begins; 

D- = subsequent delay due to one or more preemp- 
tions from service by a newly arriving serious 
case; 

N = number in system at time of routine case's 
arrival; 

= number of routine cases in system at time of 
new routine case's arrival; 

Ag = number of new serious cases arriving during 
service time of a routine case. 

a. Single Server System 

The case where there is only one server available 
is considered first. Noting that 



D 



Di + 
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for a routine case, the chance of delay can be expressed as 



P[D>0] = 1-P[D=0] 



= 1 - P [Dj_ + D 2 = 0] 

= 1 - P [D^ = 0 ,0^ = 0] 

= 1 - P[D2 = 0|D^ = 0] P[Dj_ = 0] 



Since the initial delay is exactly the delay that may 
occur in a FIFO (one priority) system, P[Dj^ = 0] is simply 
one minus the delay probability, P [T >0], as calculated in 
the previous two sections, either through use of the tables 
in Appendix A, through direct computation using Equation 5.1, 
or from the charts in Appendix B. Further, the conditional 
probability of no delay given no initial delay is 



P [D 2 = 0 1 = 0] 




fY(x) dx 



, 1/S 

^1/3 + A 



■) 



S 



k 



The equality follows from the assumption of Poisson arrivals, 
with rate A , of serious cases, and recognition of the 
integral as the Laplace Transform of the Erlang PDF f(x) . 

The final form for the delay probability is then 
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P [D > 0] 



1 

6 



= 1 






( 






-) (1 



P [T > 0] ) . 

q 



b. Two Server System 

While extension to allow multiple servers is 
beyond the scope of this thesis, a tractable result for the 
case where there are only two servers in the system is 
presented here, although the early part of the derivation 
is carried out for general c. As before. 



P [D > 0] = 1 - P [D = 0] . 



Conditioning on the total number of cases in the system at 
the time of the routine case's arrival yields 

c-1 

P[D=0] = I P[D = 0|N=i] P[N=i] (5.2) 

i=0 

The second factor is merely the steady state system size 
distribution which has been conveniently tabulated in Ref. 

4 and is contained in Appendix A. Now consider 

i 

P[D=0|N=i] = I P[D = 0|N =r,N=i]P[Nj^=rlN=i]. (5.3) 

r=0 

In this expression, the second factor under the summation 
can be recognized as the probability of r routine cases 
out of i total cases currently in service. Since the 
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probability that any one case is routine or serious is 
Ar/(Ar+As) respectively. 



P[Nr= rlN = i] 




i-r 



r 



0,1 






i = 0,1 



, . . . , 



The first factor is more difficult, however, and requires 
further conditioning on the number of serious case arrivals 
during a routine service period: 



P[D=0|N =r,N=i] = I P[D=0|a„ = j,N_ = r,N= i] (5.4) 

^ j=0 ^ ^ 

X p[Ag = j |Nr= r,N = i] . 

Since arrivals are assumed to be independent of the status 
of the system. 



Assuming the previous distributional assumptions hold, that 
is serious cases arrive in a Poisson fashion with rate A^ 
and the Erlang (k,6) distribution is appropriate to describe 
the service times. 



00 



P [Ag = j |Nr= r,N = i] = P[Ag = j]. 
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oo 



PCA =j] = / P[A^ = j |y = X] f^(x)dx 

b 0 S X 



Since 



P[Ag = j |Y = x] 



-A X (A X) ^ 
S S 

j! 



and fY(x) was given in Equation (3.1), 



P[A3 = j] 



(k+j-1) 



k j 
p 



(5.5) 



where 



= ( 



1 

B 



^ + X 
B s 



) 



and 



q = { 






) 



The resultant expression can be recognized as a Negative 
Binomial PDF although that fact will not be exploited here. 

Since this derivation is restricted to the case 
of two servers, possible values of i in formula (5.4) are 
limited to zero and one, allowing enumeration of all possi- 
ble combinations. For i = 0, r = 0 is the only possibility. 
Then 

if j = 0 , 1 
if j ^ 2 



P[D= 0|A^ = j,N_ = 0,N= 0] = 

S ^ 



r j“l 
(p ) 
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because if no serious case or only one serious case arrives 
during the routine case's service time, his service will be 
completed uninterrupted. If, however, more than one serious 
case arrives during his service time, his service will be 
uninterrupted only if the second serious case arrives after 
the completion of service of the first serious case and the 
third serious case arrives after completion of service of 
the second serious case, etc. The assumption that serious 
cases arrive in a Poisson Process at rate implies that 
the probability of a serious case's service completing prior 
to the next serious arrival is 



k 

P 



1 

6 



h 



k 



since a service time is the sum of k Exponential phases, 
each with mean duration 3. Then Equation (5.4) with i = 0 
and r = 0 becomes 



. . 1 

P [D = 0 |N_ = 0 ,N = 0] = P[A^ = 0] + I (p^)^"'^ P[A =j] (5.6) 

R S b 



Next letting i = 1 and r = 0 



j = 0 



P [D= 0 IA 3 = j ,Nj^= 0,N = 1] = 



P[S < T] (p’^) ^ ^ 



j 1 1 
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where S equals the time to completion of remaining service 
of the active serious case and T equals the time to the first 
serious case's arrival. This expression is explained as 
follows. If no serious case arrives, no interruption of ser- 
vice occurs. If, however, one or more serious cases arrive, 
the routine case's service will be uninterrupted only if 
the first serious case arrives after completion of the cur- 
rently served serious case's remaining service time and 
each new serious case arrives after the completion of 
service on the previous serious case. Therefore, Equation 
(5.4) with i = 1, r = 0 is 



P[D= 01Nj^= 0,N = 1] 



P[A_ = 0] +P[S<T] I (p^) = j] . 

^ j = l ^ 



(5.7) 



Now, letting i = 1 and r = 1, 



if j = 0 



P[D= 0|Ag = j,Nj^= 1,N= 1] = 



(ip [S > T]+P [S < T] (p^) ^ ^ if j ^ 1 



Again, no service interruptions can occur if no serious 
case arrives. If exactly one serious case arrives, service 
can continue uninterrupted only if the other routine case 
either remains in service until the serious case arrives and 
then absorbs the resultant preemption (with probability one 
half) , or if he completes his remaining service prior to the 
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serious case's arrival. All subsequent serious cases must 
then arrive after the preceeding serious case has completed 
an entire service period. Therefore, when i = 1, r = 1, 
Equation (5.4) becomes 



P[D= 0| 1,N = 1] = P[A^=0] + 

+ P [S ^T] ) 

It can be shown that 



(jP[S > T] 

I (p^)=’"^P[A =j] 
j = l s 



(5.8) 



I <5-9) 

D=1 ^ (1 - p q) 



and 



P[S <T] = 



^[1 - P ^] 



(5.10) 



The derivations are included as Appendix C. 

The means now exist to calculate the delay proba 
bility for a routine case in a two server system. First, 
however, it will be convenient to adopt the following simpli 
fying notation; 



tTq = P[N=0] 



7T- 
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P [N = 1] 



P[A =0] 
S 



k 



P 




PIA^ = jl 



k > k 
(1 - p q) 



1 



1 



j = l 



1 



P [S < T] 




Then for c = 2, Equation (5.2) yields 

P[D = 0] = P [D = 0 I N = 0]tTq + P [D = 0 I N = 1]7T^ 

Now, from Equations (5.3) and (5.6), it's easy to see that, 
for the two server case 

P[D= 0 |n = 0] = Qq 

and further that, using Equations (5.3), (5.7), and (5.8), 
P[D=0lN=l] = [Qq + Qj_Q 2 ]q + {Qq + [^d - 03 ) + Q 2 ] Qi >P 





yielding 



P [D = 0] 
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Finally then, 



P[D> 0] = 1 - P(D = 0] 



= 1 - [(Qo+Qi)^o 



I^^l (1 " §P) QiQ2^ ^1^ 



(5.11) 



c. Example 

The following example is provided to illustrate 
the computational procedure required to obtain this delay 
probability. Parameter values, though hypothetical, were 
chosen so as to provide realistic estimates of those encoun- 
tered in the course of this thesis. 

= 1.0/hr. Ag = .05/hr. 6 = .5 hrs . 




The values of the components of Equation (5.11) are provided 
be low . 

1 

p = ^ = .9756 

7Tq = . 3308 (from Appendix A) 

TT^ = .3383 (from Appendix A) 
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.9518 




k 



P 




1 



P [Ag < j ] 



.0005 




P [S < T] 



.9640 



Direct substitution yields 



P[D>0] = 1 - {[(.9518 + .0005) ](. 3303) + [.9518 + ^(.9756) 



that obtained using the same parameters but considering only 
a FIFO priority discipline. The delay probability for the 
single priority case, as taken from the charts in Appendix 
B, is only .1750, substantially less than that for a routine 
case in a two priority system. But this is to be expected, 
as it is a reflection of the opportunity for additional delay 
introduced by the possibility of bumping in the queue and 
preemption from service. 



X (.0005) + (1 - |( .9756) )( .0005) ( .9640) ](. 3383) } 



. 3629 



It is interesting to compare this result with 
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VI. CONCLUSIONS 



A. DISCUSSION 

The results presented in this thesis provide the 
machinery to conduct a sound comparative analysis of the 
merits of differing boat crew allocation schemes. Investi- 
gation of Tables III, IV, V, and VI shows quite dramatically 
the improvement in the various measures of system performance 
brought about by increasing the number of available servers. 
These results are obtainable without great computational 
capability or familiarity with the underlying mathematical 
theory. The use of one or the other of the two methods pre- 
sented is a matter of user preference. In addition, these 
techniques suggest the development of system standards by 
which proposed crew allocations could be deemed acceptable 
or unacceptable. Such a project would be a more appropriate 
task for staff members of the Search and Rescue Branch at 
Coast Guard Headquarters . 

B. RECOMMENDATIONS FOR FURTHER RESEARCH 

This thesis leaves some fairly obvious avenues for future 
work, the first being the extension to the general case of 
c servers. Perhaps, by first allowing three servers, the 
results might not only be tractable but could suggest an 
extension to allow any number of boat crews. Another area 
for further work would be to increase the number of priority 
classes. As previously discussed, boat crews are actually 
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assigned according to a continuous priority scheme, with 
each arriving case being assessed as to its severity relative 
to other active cases. While the two class case studied in 
this thesis is an improvement over a simple FIFO discipline, 
it still represents a gross simplification and could be 
improved upon. Finally, the computational techniques pre- 
sented herein raise the possibility of the creation of an 
interactive computer package designed to supply some or all 
MOE*s at the touch of a computer keyboard. Though not com- 
plicated individually, the sheer volume of computations make 
hand calculation tedious at best. 

C. PRIORITY CONSIDERATIONS 

Considering the utility of approximating the actual 
priority discipline employed by an RCC controller, and of 
applying, consequently, the methods proposed herein, the 
point of view of the SAR planner must be considered. The 
consequence, in terms of loss of life and property, of a 
failure on the part of the Coast Guard to provide a ready, 
available rescue response in the event of a serious case 
is far more grave than for that of a routine case. A fishing 
boat out of gas in a protected bay can tolerate an hour's 
delay far better than a vessel on fire or taking on water 
and in danger of sinking. Keeping this in mind, and observing 
from the data that cases of a serious nature are the rare 
exception rather than the rule, the two priority class 
assumption leads to a conservative estimate of probable 
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system performance from the point of view of the truely 
serious case. In other words, the actual likelihood that a 
distress case requiring rapid emergency response would not 
receive immediate attention will always be less than that 
predicted by the queueing techniques described in this work. 
This conservative posture, then, provides an upper bound or 
worst case prediction for the serious arrival, and therefore 
justifies its implementation until such time as continuous 
priority considerations might be properly integrated into 
the solution methods. 
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APPENDIX A 



STEADY -STATE 


TABLE I 

DISTRIBUTION AND MEAN OF NUMBER IN SYSTEM 
FOR VARIOUS M/Ej^/C SYSTEMS 



Notation : 



N 


= the number of customers in the system 

(including those being served) in steady 
state (a random variable) . 


STATE I 


= the condition of having N = I. 


P(N = I) 


= probability that N = I 


P(N < =1) 


= P(N £ I) . 


RHO 


= p, the traffic intensity. 


P( DELAY) 


= P(N^c) = P(Tg>0), the probability that 

an arriving customer cannot begin service 
immediately . 


L (GIVEN K) 


= L = E[N], the expected number of customers 
in the system. 


LQ( GIVEN K) 


= Lg = E[N] - cp , the expected number of 

customers in the queue (excluding those 
being served) . 


(LQ FOR K = 1) 


= L (k = 1) , the value of L for the corres- 

q q 

ponding M/M/c system. 


RATIO 


= L/(L^(k=D). 

q q 
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APPENDIX B 



DELAY PROBABILITY (M/M/c QUEUE) 

The following two charts provide the delay probability 
P[Tg>0] for the M/M/c queue and consequently provides an 
excellent approximation for the M/Ej^/c queue as discussed 
in Chapter V. 
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APPENDIX C 



1. Derivation of P[S^T] 

S = Excess Life of customer in service 
T = time to arrival of next serious case 



The PDF, in steady state, of the Excess Life of a non-negative 
random variable with CDF F is shown by Ross [Ref. 13] to be 
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Letting T = x and using the independence of S and T gives 
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